function prob=gampdf_approx(x1,x2,shape,scale)

prob2=gamcdf(x2,shape,scale);
prob1=gamcdf(x1,shape,scale);

prob=prob2-prob1;